arXiv:l502.00514vl [cond-mat.mes-hall] 2 Feb 2015 


The Rabi Hamiltonian in the dispersive regime 

Titus Sandu^ 

^National Institute for Research and Development in Microtechnologies-IMT, 

126A, Erou lancu Nicolae Street, 077190, Bucharest, ROMANIr^ 

(Dated: February 3, 2015) 
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I. INTRODUCTION 

The Rabi Hamiltonian was used initially to describe 
the interaction of nuclear spins with magnetic fields [l|, 
but it is also used to model the ammonia molecule 0 or 
electrons coupled with a phonon mode in a crystal lattice 
Q. More recently, beside the interaction of atoms with 
an electromagnetic field in a cavity Q the Rabi Hamilto¬ 
nian describes the interaction of superconducting qubits 
with a nanomechanical resonator [ 3 , with a transmis¬ 
sion line resonator [^, or an LC resonator Q- The Rabi 
Hamiltonian expresses the interaction of a two-level sys¬ 
tem (TLS) with a single boson mode and has the follow¬ 
ing form: 


Recently Zueco et al. [13 used the same technique 
of the unitary transformations to include the counter¬ 
rotating wave term with an improved accuracy of energy 
levels for a wider range of parameters. We build on that 
work and in Section 2 we provide not only the approxi¬ 
mate Hamiltonian and its eigenvalues but also its eigen¬ 
functions to the second order in the coupling constant. 
In section 3 we discuss comparatively the regimes of slow 
and fast TLS with respect to the oscillator and in Section 
4 we conclude the work. 


II. THE APPROXIMATE HAMILTONIAN IN 
THE DISPERSIVE REGIME 


H = uta'a -b Hctz + A(a^ -|- a){a- -b tT+), (1) 

where H = 1, ax, ay, and CTz are spin 1/2 matrices, 
a± = CTa; ± iay, and a are the creation and annihi¬ 
lation operators of the quantum oscillator. 

The apparent simplicity of © led to the conjecture 
about its solvability Q that only recently has been fully 
proved Despite the fact that the Rabi Hamiltonian 
is completely solvable, its general solutions cannot be 
easily grasped in simple terms. In quantum optics with 
atoms in a cavity also known as cavity quantum elec¬ 
trodynamics (cavity QED) one encounters a regime of 
quasi-resonance (w ~ H) and a coupling strength X/to 
between 10“^ and 10“^. This regime is very well de¬ 
scribed by the Jaynes-Cummings model which is simpler 
and solvable [l^. The Jaynes-Cummings model uses the 
rotating wave approximation (RWA) in ([T]) , such that the 
counter-rotating wave term a -b aa- is cornpletely ig¬ 
nored 0. Moreover, in the circuit QED [10, El , where 
the TLS is associated with a qubit, one encounters a dis¬ 
persive regime (i. e., large detuning |a; —D| ^ w-bD) with 
a ultra-strong coupling A/ a; ^ 0.1, where the Jaynes- 
Cummings model fails |12l - H4| . A more elaborate ap- 
proximation based on the polaron transformation [^, [l^ 
is the generalized RWA, which works well in the region of 
zero and large detuning and ultra-strong coupling (T^ . 


‘Electronic address: titus.sandu@imt.ro 


Quite often a general Hamiltonian H can be decom¬ 
posed as H — Hq + X W, where Hq is diagonalizable and 
A is a small real parameter. The ultimate goal of the uni¬ 
tary transformation approach is to find an anti-hermitian 
operator S such that the unitary transformation U = 
diagonalizes H. However, the complete diagonalization is 
not always possible but it can be done successively begin¬ 
ning with the first order of A. The requirement to elimi¬ 
nate XW in the first order leads to the form of S that is 

0 

given by S' = —i A lim J e^*Wi (t) dt, where Wj (t) = 

^^°+-oo 

choose Hq = oja^a + Aaz and 
XW = X {a+ a^) {a- + CT-i-) we can calculate S as being 


S = - 


2iu! 

A2-a;2 





2iA 

A2 -w2 


(a -b a^) 


( 2 ) 

The use of © in the transformation Ui = changes 
the Rabi Hamiltonian into 


- a. . 2A^A a a., 2 A^W / a ■? \ 

H — uja^a+Aaz + -^ -y (a + o') + ) ■ 

( 3 ) 

We notice that the small parameter is A = }? j { A? — 
and if A <C 1 we can retain just the first three terms of 
@. Then we expand (a -b a^) , make some arrangements 
and remove completely the part of Hamiltonian propor¬ 
tional with . This is possible with the unitary 

transformation U 2 = \ where /3 = 


/3+ 0 
0 /3_ 
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and tanh{2l3±) = . Thus, the approximate form of (jS)) is 

I 


T, + r~n It I V + SAwAcr^ — w ,, . t 

Hi = a)a\JuP- + SAwAtr^ H-^-f (1 + 2A)Acr2 + Aw. (4) 


To be fully consistent with the approximations made so 
far we expand the radical in Eq. (j3]) and keep just the first 
two terms in the series. Thus, the oscillator frequencies 
are uj± = w ± 2AA and the Hamiltonian (|3]) turns into 
the Hamiltonian of the dispersive regime and ultra-strong 
coupling 

Happrox — ^0 w ^ dAjAfj^ “t“ Aw. (5) 

It is easy to see that the eigenvalues of are 

E±n = nui± ± ——A -f Aw. (6) 

I 


We also calculate the eigenvectors of ([5]) in the original 
Schrodinger picture. In order to do so we need to evaluate 
the action of on any vector in the Hilbert space of 
the problem. The action of is rather cumbersome 
but a meaningful expression can be obtained with the 
Zassenhaus formula (l9l |. which, for any operator A and 

B reads — q\a^xb^-^[a,b] . Thus, the 

eigenvectors of are 


|«-±„ >= e 


■ 2XA 

“ “'e 






(7) 


In Eqs. ® and 0 n is a natural number, |n(w±) > 
are the nth eigenvectors of the quantum oscillator with 
frequencies w±, and s± are the eigenvectors of Uz, i. e., 
azS± = (±I/2)s±. 


III. DISCUSSIONS 

The energy levels given by Eq. ® have been compared 
with the exact eigenvalues of o in the paper of Zueco 
et al. [13, where it has been proved a very good match 
between dS]) and the eigenvalues of ([T]) for a wide range of 
parameters in the dispersive regime. On the other hand, 
the accuracy of eigenfunctions has been checked in Ref. 

by calculating their fidelity with respect to the exact 
numerical calculated eigenfunctions. The comparison of 
the eigenfunctions have been performed in two limiting 
cases of the dispersive regime: slow TLS and fast oscilla¬ 
tor (A "C w) and slow oscillator and fast TLS (A u). 

In the slow TLS and fast oscillator regime (A <C w) 


Agarwal et aL[2^ have found that the eigenvectors of 
Hamiltonian ([S]) in Schrodinger picture have a poor fi¬ 
delity with respect to the exact eigenvectors. Moreover, 
their calculations indicated that the fidelity of the eigen¬ 
vectors provided by the adiabatic approximation |15l . Il6l | 
is excellent. We compare our results given by ([7]) with the 
adiabatic eigenvectors. We recall that the adiabatic ap¬ 
proximation is obtained using the polaron transformation 
neglecting the spin non-diagonal term 
[l5l |. Hence, the adiabatic Hamiltonian and its eigenvec¬ 
tors are, respectively. 


Hadiab = wa'^a -I- /\azCos{ 


2\a — a^^ A 


( 8 ) 


\^±u.ad^ab >= “^)'"-|n(w) > s±. (9) 

When A <C w both A/w and A/w are small, therefore d?]) 
becomes 


J 


I- 


( 10 ) 


In Eq. (nni the part that contains the squeeze operator deviates from the unit operator by an amount of the same 
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order as the overlap < n(a;)|n(a;±) > differs from unity, 
i. e., thus their deviation from 1 can be safely 

discarded. It is easy now to evaluate the fidelity between 
dH) and (Uni) as / = I < 'i’±n.adiab\'i’±n > P = 
which is very close to 1. It would imply that the fid elity 
of (O is close to one and in stark contrast with Ref. . 
One possible explanation of this discrepancy might be 
the way is calculated to obtain the eigenfunctions. 
We adopted an exponential approximation based on the 
Zassenhaus formula, but if one adopts a series expansion 
of the fidelity would be / = 1 — 0{X/uj), which is 
sensibly below unity and consistent with the conclusions 
of Agarwal et al. 

In the region A:^ uj (fast TLS and slow oscillator) one 
can also invoke an adiabatic approximation as long as the 
TLS and the oscillator are on two different energy scales. 
This situation is more often encountered in molecules and 
crystals 0. The new adiabatic regime can be obtained 
by a unitary transformation that also works well for the 
regime of deep ultra-strong coupling when X > uj [2l|. 
This adiabatic approximation generates an adiabatic po¬ 
tential with two sheets. The lower sheet has two minima 
if A > '/ujA and only one minimum if A < \/ujA [^ . 
When A > \/uiA an approximate solution based on the 
displaced oscillators generated by the the adiabatic po¬ 
tential is given in Ref. [1^. However, the case of the 
adiabatic potential with just one minimum can be de- 
scribed_W the Hamiltonian ([S]) of which both the eigen¬ 
values [l3| and the eigenvectors are well reproduced 
by Eqs. dSD-©. We note here that in the regime A ^ ui 


a series expansion of gives rather similar results with 
Eq. ©. Moreover, one can easily check that the curva¬ 
tures of the two adiabatic sheets in the origins are just 
and respectively. The eigenvector © exhibits 
also a certain degree of squeezing. In f act, the largest 
squeezing is encountered at A = Ac = \/wA where 
the system undergoes a sharp transition [2^. 


IV. CONCLUSIONS 

In this work we studied the Rabi Hamiltonian in the 
dispersive regime and ultra-strong coupling. It was used 
a unitary transformation that takes into account terms 
beyond the rotating wave approximation. We are able to 
build the approximate Hamiltonian with its energy levels 
and its eigenfunctions. It turns out that an exponential 
approximation of the eigenfunctions is better suited than 
the approximation made by series expansion. We also 
compare and discuss this approximation with respect to 
other approximations in the regimes of fast and slow two- 
level system. 
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